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We argue that dynamics of gapless Fractional Quantum Hall Edge states is essentially non-linear and 
that it features fractionally quantized solitons propagating along the edge. Observation of solitons 
would be a direct evidence of fractional charges. We show that the no n- linear dynamics of the 
Laughlin's FQH state is governed by the quantum Benjamin-Ono equation. Non-linear dynamics 
of gapless edge states is determined by the nature of gapped modes in the bulk of FQH liquid. The 
dispersion of edge modes is traced to the double boundary layer of FQH states. 
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I. Introduction and Results. In a Fractional Quantum 
Hall state electrons collectively constitute an incompress- 
ible liquid almost free from dissipation. Excitations in 
this liquid are gapped by an energy A„ determined by 
the Coulomb interaction 1 . The gap is large compared 
with the temperature (but small compared with the cy- 
clotron energy A v <C Huj c ). The only low energy current 
carrying states in a FQH liquid are edge states localized 
on the boundary [2] . Edge states provide a valuable tool 
to probe FQH states. 

In most FQH states the excitation gap is large com- 
pared to the energy of long wave edge states, and for that 
reasons is commonly neglected. The standard approach 
to the theory of Edge states starts from the Chern-Simons 
action in the bulk [3]. This action has no scale. It ne- 
glects gapped bulk modes but focuses on braiding proper- 
ties of FQH states. The Hamiltonian of this theory differs 
from zero only by a confining potential. In the case of a 
"sharp" boundary it is a topological field theory where 
all states are in a correspondence with boundary states. 
They are chiral bosons with the current algebra 



2u. 



[f(x) > f H (x')] = -V x 6(x-x% 



governed by a linear wave equation 



Here 



/ - c V x f = 0. 
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is the Hilbert transform of f(x), Ib is the magnetic length 
2n£ B = &o/B, c = h~ 1 i 2 B \V y U\ is velocity of sound 
determined by the slope of the confining potential U (y) . 
The unperturbed droplet occupies the half plane y < 0, x 
is a coordinate along the boundary. 

This theory assumes that a FQH state does not change 
toward the boundary and neglects electric polarization 
caused by edge waves. This happens if the curvature of 
the potential is small compared to the gap <g\7j[/ « A„, 
but a slope is larger than electric field £ B V y U 3> e 2 . 



Also it assumes that boundary waves are long and their 
amplitudes are small compared with scales set by the 
magnetic length. We accept these conditions. 

Physics missed by this otherwise successful theory can 
be seen in the following setting. Let us suddenly perturb 
the edge by a classical instrument, say RF-source (whose 
spatial extent is larger than magnetic length Ib) and then 
release the system. A smooth semiclassical density profile 
fo(x) will occur. How does the density profile propagate 
along the edge? The wave equation (2) suggests that 
the initial profile simply translates as f(x,t) = fo(x — 
Cot) without changes. This may be true shortly after the 
perturbation, but at time of order of h/A u the profile is 
expected to change. This time is not too long. In typical 
v = 1/3 samples it is about lOps pQ. At that time the 
wave equation fails. 

Another apparent problem of the wave equation is that 
it does not discriminate between fractional and integer 
QHE where a mobility gap is due to a disorder. Do these 
very different cases share the same dynamics? 

Corrections to linear waves (|2| come from few sources: 
the curvature of the confining potential at the edge, mix- 
ing with higher Landau levels, disorder, and, more in- 
terestingly, the interaction between the gapless edge and 
gapped bulk modes. The latter is the most relevant for 
FQHE. It is the subject of this letter. 

We show that the wave equation receives important 
corrections of two sorts: a nonlinear term and a dis- 
persion term. Both are proportional to the scale k ~ 
A v £ B /h of gapped excitations omitted in the Chern- 
Simons action. We will show that they are relevant for 
the non-equilibrium physics of FQHE at the edge. 

In this paper we focus only on the simplest Laughlin's 
states with a single branch of excitations assuming that 
the state remains intact all the way to the edge. In this 
case the quantum equation for edge modes reads 

/- c V x f -kvJ^-t)- V x f H ^j = 0, r, = ^ (4) 

Here rj = is the dipole moment of the boundary dou 



ble layer explained below (see (28) and Fig. [2]). 

The new terms in brackets cause new phenomena. One 
is fractionally charged solitons on the edge. 
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FIG. 1: Boundary waves: the boundary layer is highlighted 
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The chiral boson in this equation has a transparent 
interpretation. Acting on coherent state with an elec- 
tronic density p(x, y) the chiral boson means a boundary 

density f(x) = Jq p(x, y)dy, where y(x) is a bound- 
ary displacement counted from an unperturbed bound- 
ary y — 0. The bosons act in the chiral Fock space, where 
\hnf H + vSI x tt f—0. Here 717 = ~^J] ^ s a momentum of 
/. The constraint yields the algebra ([TJ). 

The non-linear equation Q previously appeared in two 
domains of physics. Its classical version has been derived 
by Benjamin [5] in 1967 for inner waves in a deep strat- 
ified incompressible fluid with a rapidly changing den- 
sity or shear [7]. It is called the Benjamin-Ono equation. 
The quantum Benjamin-Ono equation (qBO) identical to 
( 2p l describes the chiral sector of Calogero model [H[S]. 

Both relations are not accidental: a FQH state is an 
incompressible quantum fluid with rapidly changing den- 
sity and shear at the boundary as shown in Fig. [2j Both 
FQH and Calogero liquids feature excitations with a frac- 
tional charge. One can treat the result of this paper as a 
proof of a long anticipated equivalence between Calogero 
liquid and FQH Edge states. 

The qBO has an intrinsic relation to a boundary Con- 
formal Field Theory with the central charge c = 1 — 
Qu^ 1 {u — l) 2 situated in the exterior of the droplet. Con- 
formal symmetry emerges with respect to deformations 
of the boundary of a FQH droplet. Also, qBO has a natu- 
ral extension to non-Abelian FQH states. We will discuss 
these aspects elsewhere. In this paper we mention a few 
major features of the qBO equation jU [5] and focus on 
its derivation. 

2. Solitons in Quantum Benjamin-Ono Equation. The 
qBO is a Hamiltonian equation. The dispersion term 
V^xf H nas the dimension of viscosity but contrary to 
real viscosity does not produce dissipation. It gives a 
non-analytic correction to the dispersion of linear waves 



H(uj(k) — cok) = K7]k\k\. 



(5) 



A noticeable feature of the qBO is that the "dissipation- 
less viscosity" rj is expressed solely in terms of the filling 
fraction v and, in this sense, is universal. Conversely in 
classical liquids [BJ [7] the coefficient r\ depends on the 
equation of state. It is not quantized. Also, n vanishes 
for the IQH. It is inherent to FQHE. It is similar in origin 
to the "Hall viscosity" [9], but is not directly related. 

Another noticeable feature is that qBO being a non- 
local equation, nevertheless is integrable [8] . In integrable 
equations a competition between non-linear and disper- 
sion terms yields solitons. In the qBO there are two 
branches of solitons: one is ultrasonic, another is sub- 



sonic. Quite remarkably, both carry quantized electron 
charges [4]. An ultrasonic soliton carries an integer of 
electron charge q = +e. It is a bump on the edge - a 
coherent state of an electron. A subsonic soliton is a 
coherent state of holes - a dent on the edge. It carries in- 
teger of a fractional charge q = — ve of an opposite sign. 
Shapes of the elementary solitons are especially simple: 



f q (x + c t) = 



TT a 2 {x — Vqt) 2 + 1 



q = 1, -v. (6) 



Velocity of a soliton (relative to the sound) is v q = qnna is 
inversely proportional to the magnetic field. It is propor- 
tional to its amplitude a and its charge q. The amplitude 
a > is arbitrary, but the charge is quantized. 

Notice that in the classical limit v 0, the fractional 
("hole") branch disappears, but the quantized "parti- 
cle" branch remains. BO is the unique classical equation 
which has quantized solitons. Charge carried by solitons 
of any other soliton equations can be arbitrary. 

When a dent in the boundary density fo(x) is created 
by an RF source, corrections to the wave equation be- 
come important at time t c ~ (k|V/o|) , which is only 
by a factor {Info) ~ 1 larger than the typical time H/ A„ ~ 
lOps. At that time a wave packet collapses through a 
shock wave to oscillatory features which further separate 
to a stable fractionally charged soliton train - a sequence 
of pulses with fractionally — ve charges 4J. This suggests 
an appealing prospect for a direct observation of frac- 
tional charges: apply an RF-source to the edge, observe 
a sequence of fractionally charged pulses. 

3. Phenomenological Hamiltonian is the starting point 
of the analysis. But first we must identify a space where 
the Hamiltonian acts. This space is the result of projec- 
tion enforced by the condition Hlj c 3> A„. It is the set of 
states obtained by a deformation of the Laughlin ground 
state ipo by holomorphic polynomials. In a radial gauge 
suitable for a central-symmetric confining potential co- 
herent states of this space are 



(7) 



where ^ = Yli>ji z i ~ z j) an d Z v is a normalization. 

A complex potential V{z) is analytic at infinity and 
such that Atto — —AV is real (Imcr = 0). A meaning of a 
is a density of " holes" . From hydrodynamics perspectives 
they are vortices, cf. ( |14[ ). Also a set of permissible 
operators is spanned by a product of holomorphic and 
anti-holomorphic operators. 

We denoted the averages of symmetric operators in a 
given Estate as (0) v = J V v Oip v Ui^n. 

These states have been studied in [i"2] . We mention 
an important sum rule that follows from the value of the 
dilatation operator ' v + ' ^iV)v = 2 



N ^ 
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NP - ^{n ■ ViRevO =l-f. (8) 



3 



We construct the Hamiltonian based on a few defining 
properties: (i) Lauglin's w.f. is the ground state; (ii) the 
Hamiltonian density is an ordered product of holomor- 
phic and anti-holomorphic operators; (iii) long waves of 
a FQH liquid are Galilean invariant; (iv) on closed man- 
ifolds (no boundaries) all states are gapped [13]. 

Under these assumptions the zero Hamiltonian of the 
Chern-Simons theory is replaced by the Bernoulli energy 



H 



lv ~+ ,2 

— v'pvd r, 



(9) 



where m w = 



/<-/ 



— is an effective mass obtained 



from the value of a gap, p{r) = J2i S( r ~ r i) is the density, 
and v = i> x — iv y is the complex velocity. Velocity of a 
particle is 
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P=\- (10) 



Here A = A x —iA y is an external e.m. potential B = V x 
A, Zi = Xi + iyi is the complex coordinate of a composite 
particle, 2d z = d x — id y . We occasionally use /3 — \ jv. 

A subtle point of this Hamiltonian is the definition 
of the velocity operator v. It differs from the velocity 
of individual electrons but corresponds to the velocity 
of "guiding centers" or "composite particles" - electrons 
with an attached flux converting them to bosons. The 
velocity of guiding centers changes slowly in long-wave 
excited states. This is the velocity entering hydrodynam- 
ics. We emphasize that the "mass velocity" (10 1 differs 



from the velocity obtained from electric currents carried 
either by particles or by guiding centers |11) . 

The Hamiltonian ([9]) can be treated as a quantized ver- 
sion of the " effective" Hamiltonians proposed and studied 
inRefs. [TO]. 

4- Quantum hydrodynamics of the FQHE liquid de- 
scribes dynamics of velocity and density fields when the 
number of particles is large. The fields are defined as op- 
erators acting on averages (O) v . In particular, the veloc- 



ity field is defined as p(r)v(r)(0) v — 8(1 



)v 3 -0), 



where p(r) — (/5(r)) 1 
reads 



In this representation the velocity 



ra„v = 2d z (n p - iV), tt„ = -ih— . (11) 

bp 

In the restricted space (|7|) matrix elements of ir p are imag- 
inary, AV is real, hence the liquid is incompressible 



V • v = 0. 



(12) 



It is customary descryibe the incompressible flow in terms 
of the stream function 



v = V x <I>, = yk p 



ReV. 



(13) 



Eq. ( 13 ) gives an interpretation to the deformation po- 



function. Diagonal part of vorticity and energy operators 
f9fll| ) are 



2 m,, 



|W|>d r, m„<V x v) v = -AV (14) 



Hydrodynamics of the FQHE fluid can be also seen as 
dynamics of the function V. 

5. Holomorphic fields, potential incompressible flow 
and Edge states. In a system without boundaries all 
modes are gapped. If there is a boundary, gapless edge 
states emerge. In order to focus on edge states it is suf- 
ficient to consider only a potential flow where a stream 
function is harmonic 



0, A* = 0. 



(15) 



Potential flow corresponds to deformations of the w.f. ^ 
by analytic functions inside the domain occupied by the 
liquid. All singularities of V are outside of the domain 
(analytic functions do not exist on closed manifolds, so as 
gapless modes) . In this case a FQH flow is potential and 
It has been shown in [T2] that a holomorphic deformation 
of Laughlin's state changes only the shape of the droplet, 
leaving the density and the area unchanged. In the lead- 
ing 1 /N order, the bulk density is uniform p = i/B/Qq. In 
a central-symmetric potential the droplet in the ground 
state is a disk with a radius R = \J N/np. 

Incompressible potential flow with a free boundary 
and a constant density is a standard subject in classi- 
cal hydrodynamics [14] . Its extension to the quantum 
case is straightforward. Use the Green formula and the 
Cauchy-Riemann condition to express the Hamiltionian 
only through the boundary value of the fluid potential 

H=^<f VdnVds, (H) v = -B- I VdV. (16) 
2 J 2im v J 

The bulk Hamiltonian vanishes. 

The remaining task is to connect the potential V, or 
the velocity to the boundary elevation yix). Then the 
only governing equation is the kinematic boundary con- 
dition [14] 

V + v x V x y + v y = 0. (17) 
Later in this paper we show that 

v x = c - Kpy(x), v y = Kn- y^ x (18) 



where 



are velocities of the inner layer tangen- 



tential V. Its real part is the diagonal part of the stream 



tial and normal to the unperturbed boundary. Combing 
( 17|18 | and using f(x) = py{x) we obtain the qBO Q. 

6. Chiral constraint is a relation (JTsJ) between the ve- 
locity and the shape of the boundary of the droplet, hence 
a relation between V and y(x) in a state ([7]). This relation 
has been studied in [T^]- Here we obtain the constraint 
invoking Dyson's arguments [15] and refine the results 
of [12]. Dyson's arguments are somewhat heuristic but 
transparent and short. They are somewhat heuristic but 
transparent and short. 



4 



Let us express an expectation value (0) v as a path 
integral over the density field 

(0) v = [ 0[p]e-^Vp. 



The chiral constraint (aka loop equation |12j ) is the saddle 
point condition ensured by a large number of particles 
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f3-F v [p])0[p] 



0. 



(19) 



The functional F v [p] can be treated as the free energy of 
2D-Coulomb plasma. It consists of energy and entropy. 



0F V =loghM 



plogpd r. 



(20) 



The entropy is the Jacobian of passing from integration 
over particle coordinates to a path integral over the den- 
sity field. 

In order to find the energy of the plasma we write 



log | 



= E I >g| 



where i(r) is the mean distance between particles. Ex- 
clusion of " self-interaction" allows to replace sums by in- 



tegrals: Epilog 



4v( r i) ~ log^(rj), where 



(p(r) 



2 / log \r — r'\pdr 



(21) 



is the Coulomb potential of the plasma. This gives 
PF V = f[^(tp-tp)+^\ogl + \ogp-^ReV}pd 2 r, (22) 

where tp — — 2 ^ log \r — r'\pd 2 r' = irpr 2 is the potential 
of a neutralizing uniform background charge p. 

A subtle point of this approach is the value of the mean 
distance between particles. It is different in the bulk and 
close to the boundary. In the bulk the mean distance is 
isotropic and £ ~ 1/ ^fp. The short distance term and the 



entropy term in (22 1 sum up to (§• — l)logp 



Summing up we write the condition ( 19 1 as 



bulk: 2ttk x ^ = ip - ip - (1 - 2v) log p. 



(23) 



Close to the boundary, i entering (22 1 is the distance 



in the direction normal to the boundary. Since the mean 
distance along the boundary is constant ~ £b we have 
IbI ~ In this case the short distance term and the 
entropy term sum up to (/? — 1) log p. 

Unfortunately, Dyson's arguments miss some exponen- 
tial corrections which become important at the bound- 
ary. Notice that in the case of the IQHE when j3 = 1 the 
term (/3 — 1) logp vanishes. The Dyson's arguments give 
2itk~ 1 '$> = tp — <p . This equation treats the density as a 
step-function. Instead, the exact density of the IQHE is 



Pi 



2 p_erfc( 7fc : 



j(e(- tf ) 



>2lT 




FIG. 2: Boundary double layer of v = | state: Left (in de- 
creasing order): density p, density pi at v = 1 (scaled), the 
dipole moment —Tj(r) = —f r y'(p — pi)dy' . The value of 
|oo = ^Er ~ 0.05307 is clearly seen. Right: p— pi illustrates 



Eq.(30l [19] ■ The inner boundary lays in the bulk before the 
overshoot. 



where (— y) is the distance to the boundary. 

This failure can be "repaired" by replacing the plasma 
Coulomb potential tp = irpr 2 of a neutralizing uniform 



charge p in (22) by the potential of the charge pi: tp/ = 
-2/log|r--r'|p / (r')d 2 r' « npr 2 (l + 0(e-^ '^l)). Then 

boundary: 27tr~ 1, I' = tp — tpj — 2(1 — v) log p. (24) 

This ad hoc procedure reflects a disreteness of particles. 

So far we did not assume that the flow is potential. If 
the flow is potential the Laplace operator nulls the l.h.s. 
of (23). We obtain a Liouville-type equation 



boundary: p — pi + 77A log p = 0. 



(25) 



We do not have a satisfactory mathematical justification 
of this equation, but it treats correctly exponential cor- 
rections becoming important at the boundary. Also, nu- 
merical solution of the Liouville equation ( [25] ) appears to 
be undistinguishable from a numerical ah initio simula- 
tion of the Laughlin's wave function ([7|) [18] [19] . 

Confining potential is easy to incorporate into ( 24 ) . It 



shifts the current pv by the Hall current x U . Taking 
the curl of ( 24 ) and multiplying by p we obtain [16] 



— (pv-^VxlT) = n^xUp-tpA- 2(l-i/)Vxp (26) 
K h 



We assume that l 2 B V 2 y U < A„ 



Then the effect of the 
confining potential at the inner boundary is reduced to 
a shift of the velocity along the boundary v x — > v x — cq. 
Eqs. ( [23] 26] ) can be checked against the sum rule p[). 

Using ill) and an exact formula J r 2 pid 2 r = ^-(JV+l). 
we can write the sum rule as 



1 



d 2 r[f x Gov)] — 



ttR 2 



r 2 (p- Pl )d 2 r = 477. (27) 



7. Boundary double layer and its dipole moment Eq.(25) 
has a profound consequence. It shows that the density 
behavior is singular: on approaching the boundary the 
density oscillates and shoots up before falling down. The 
overshoot has been observed numerically in |17j and also 
discussed in [T^] . A detailed structure of the overshoot is 
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complicated and not yet satisfactory understood. How- 
ever, we know that in a rough but sufficient approxi- 
mation it is a double layer as is in Fig. [2] As we see 
in a moment, only the Coulomb potential of the layer - 
the boundary dipole moment (per unit length) enters the 
edge dynamics. 

The value of the dipole moment follows from the sum 
rule ( 8|27 l evaluated sufficiently close to the boundary 
where normal velocity in (27 1 can be neglected. It is 



1 



4?r 



(28) 



A detailed structure of the boundary layer follows from 
(25 1. Iterations of this equation allow us to conjecture 



the asymptotic expansion for the density from the bulk 
side. In units of Ib it reads 

-£ 2 /2 

V p = p(y+ e —^{2^r,i + 0{r 1 )), £<0. (29) 



At Ib — > the leading term of ( 29 1 is a double layer 



p(y) ~ pi(v) + v$ l (y)- 



(30) 



The following arguments are based on this major formula. 

8. Transformation of velocities Now we are in a po- 
sition to compute the velocity in terms of the bound- 
ary elevation y{x) and to complete the governing equa- 
tion (17 1. As in (17) we assume the density moves to- 
gether with the boundary p(x,y) — po(y — y(x)), where 
po (y) is the density of the ground state and compute the 
boundary value of the Coulomb potential ip ( [21] ) enter- 
ing the chiral constraint ( 24 ) . This exercise is equivalent 



to the Hadamard formula for a variation of solution of a 
boundary value problem upon a variation of the bound- 
ary. Computing / log \r — r'\po(y' — y(x'))dx'dy' , we shift 



the variable y' y' + y{x'), subtract pj and expand in 
y(x). We obtain ip(x,y) - Lp (y - y(x)) w 



dx'(y(x')-y(x)) 



y-y 

r - — 2 



[po{y')-pi\ d y'-^pyy{ x )- 



The integral over y' is localized inside the boundary layer. 
If we choose y to be on the inner boundary of the layer 
the range of \y — y'\ ~ Ib- We can safely replace \r — r'\~ 2 
in the integral by [x — x)~ 2 and obtain the transforma- 
tion law of the stream function and the velocity under 
a displacement of the boundary by y(x) (valid for any 
shape of the boundary) 

#(sb, y) = # (v - y{x)) -n(py v( x ) + (31) 



Vox 



Kpy(x), 



H 



'Jy = VOy + IVr) ■ Vxx- 



(32) 



For a flat boundary vq x — cq, vq v = 0. 



A relation between V and y(x) follows from (32) 
p f y(x')dx' rj f y(x')dx' 



4tt 



2tt 



{z-x'f 



(33) 



The Benjamin-Ono equation Q and the current al- 
gebra ([I]) follow from the transformation law ( 32 ) , and 
the value of the dipole moment (28). Alternatively the 



Benjamin-Ono equation can be also obtained from the 
boundary Hamiltonian ( 16 1 expressing it through / with 
the help of pi. 
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